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RHIC experiments have recently measured the azimuthal correlation function of forward di- 
hadrons. The data show a disappearance of the away-side peak in central dn-Au collisions, com- 
pared to pH-p collisions, as was predicted by saturation physics. Indeed, we argue that this effect, 
absent at mid-rapidity, is a consequence of the small-x evolution into the saturation regime of 
the Gold nucleus wave function. We show that the data are well described in the Color Glass 
Condensate framework. 
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Hard processes in hadronic collisions, which resolve the partonic structure of hadrons, are 
well described by the leading-twist approximation of QCD. In this weak-coupling regime, partons 
in the hadronic wave function scatter independently, this is the essence of coUinear factorization. 
However, since the parton densities grow with decreasing energy fraction x, the hadronic wave 
function also features a low-x part where the parton density has become large and partons scatter 
coherently, invalidating the leading-twist approximation. This weak-coupling regime, where non- 
linearities are important, is called saturation, and it can be probed at high-energies since increasing 
the energy of a collision allows to probe lower-energy partons. 

In hadron-hadron collisions, by contrast with deep inelastic scattering, hard processes are sin- 
gled out by requiring that one or more particles are produced with a large transverse momentum, 
much bigger than Aqcd- When in addition the particles are produced at forward rapidities, such 
processes are sensitive only to high-momentum partons inside one of the colliding (dilute) hadron, 
whose QCD dynamics is well understood, while mainly small-momentum (small-^c) partons inside 
the other dense hadron contribute to the scattering. Replacing that hadron by a large nucleus further 
enhances the gluon density, and the possibility to reach the saturation regime. 

In the case of single-inclusive hadron production, the suppression of particle production at for- 
ward rapidities in d-i-Au collisions compared to p-i-p collisions, experimentally observed at RHIC 
[|l|], constitutes one of the most compelling indications for the presence of non-linear QCD evo- 
lution effects in presently available data. The Color Glass Condensate (CGC) provides a robust 
theoretical framework to describe the small-x degrees of freedom of hadronic/nuclear wave func- 
tions. The good description of, among other observables, forward hadron production in d-i-Au 
collisions at RHIC indeed lends support to the idea that saturation effects may be a relevant 
dynamical ingredient at present energies. 

However, alternative explanations of the suppressed forward hadron yield in d-i-Au collisions 
were proposed |^, suggesting that one is not yet sensitive to the saturation of the nuclear gluon 
density at RHIC energies, but that the suppression is rather due to partonic energy loss through the 
nuclear matter, neglected in CGC calculations. In spite of the fact that saturation-based approaches 
were the only ones to correctly predict this phenomenon, the existence of alternative scenarios 
calls for the study of more complex observables. In the light of recent preliminary data in d-i-Au 
collisions at RHIC, showing the production of forward mono-jets calculating double-inclusive 
forward hadron production in both frameworks could help pin down which is the correct picture. 
In this work, we show that the CGC calculation predicts correctly the azimuthal de-correlation 
of forward di-hadrons in d-i-Au collisions compared to p-i-p collisions [Q], thus providing further 
support for the presence of saturation effects in present data. 

1. Formulation 

In the case of double-inclusive hadron production, denoting p\±, p2± and yi,y2 the transverse 
momenta and rapidities of the final state particles, the partons that can contribute to the cross 
section have a fraction of longitudinal momentum bounded from below, by Xp (for partons from the 
deuteron wave function) and xa (for partons from the nucleus wave function), which are given by 

Xp =x\ +X2 , xa=x\ e^^-^'' +X2 e^'^^^ , with x, = e^' . (1.1) 
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The kinematic range for forward particle detection at RHIC is such that, with y/s = 200 GeV, 
Xp ~0.4 and xa ~ 10^^. Therefore the dominant partonic subprocess is initiated by valence quarks 
in the deuteron and, at lowest order in a^, the dAu^h\h2X cross-section is obtained from the 
qA — )• qgX cross-section, the valence quark density in the deuteron f^i^, and the appropriate hadron 
fragmentation functions D/j/^ and Di^jg. 



A, Jx2 J-^+-^ 



dN'i^^''^''(xPA,-)DH,/,{zu^)Dh2/giz2,^)+ 

\ Z\ Z2 J 



Z2 Zl 



Here we will use the CTEQ6 NLO quark distributions and the KKP NLO fragmentation functions. 
The factorization and fragmentation scales are both chosen equal to the transverse momentum of 
the leading hadron, which we choose to denote hadron 1, jj. = \pi±\. Note that we have assumed 
that the two final-state hadrons come from partons which have fragmented independently, therefore 
formula (L2) cannot be used when = (j2 — Ji)^ + (A^)^ is too small, where is the differ- 
ence between the hadrons azimuthal angles. Computing the cross section for small R requires the 
introduction of poorly-known di-hadron fragmentation functions, we shall not do it, because as we 
shall see the non-linear QCD effects we are interested in manifest themselves around AO = n. 

As usual, due to parton fragmentation, the values of x's probed are generically higher than 
Xp and Xa defined in ( |1.1[ ). For the proton, one has Xp < x < I, and if Xp would be smaller (this 
will be the case at the LHC), the gluon initiated processes gA — qqX and gA — ggX should also 
be included in (|L^), they have been computed recently [^] . For the nucleus, we shall see that 
the parton momentum fraction varies between xa and e^^-'^ +e^^^'^. Therefore with large enough 
rapidities, only the small-x part of the nuclear wave function is relevant when calculating the qA — )• 
qgX cross section, and that cross section cannot be factorized further: dN'^^^'^^^ / fgjA^dN^'^^^'^. 
Indeed, when probing the saturation regime, dN'^^^'^^^ is expected to be a non-linear function of 
the nuclear gluon distribution, which is itself, through evolution, a non-linear function of the gluon 
distribution at higher x. 

Using the CGC approach to describe the small-x part of the nucleus wave function, the qA — 
qgX cross section was calculated in ||7|, ^]. It was found that the nucleus cannot be described 
by only the single-gluon distribution, a direct consequence of the fact that small-x gluons in the 
nuclear wave function behave coherently, and not individually. The qA — )• qgX cross section is 
instead expressed in terms of correlators of Wilson lines (which account for multiple scatterings), 
with up to a six-point correlator averaged over the CGC wave function, while the gluon distribution 
is the Fourier transform of only a two-point (dipole) correlator =yK: 

where r denotes the dipole transverse size. F{x,k±) is actually called the unintegrated gluon distri- 
bution, due to the fact that it is k± dependent, a feature known to be necessary to describe small-x 
partons, even in the linear regime. 

At the moment, it is not known how to practically evaluate the six-point function. In [||], 
an approximation was made which allows to express the higher-point correlators in terms of the 
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two-point function (1.3). This is done assuming a Gaussian distribution of the color sources, with 
a non-local variance. The resulting cross section for the inclusive production of the quark-gluon 
system in the scattering of a quark with momentum xP^ off the nucleus A reads ^ : 



2 



, (1-4) 



where q and k are the momenta the quark and gluon respectively, and with A = k± + q± and 
z = k^ /xP^ . In this formula, xa denotes the longitudinal momentum fraction of the gluon in 
the nucleus, and xa = xi e^^^'/zi +-^2 e^^^^/z2 > xa when the cross section (O) is plugged into 



formula (1.2). 



The second line of formula (1.4) features the so-called A;/- -factorization breaking term, with 



I^p{z,kr,x) = J d^q^\l/^p{z,qi_)F{x,k^-q^) , (1.5) 

and where Yap is the well-known amplitude for q^qg splitting (A , a and j3 are polarization 
and helicity indices). While no additional information than the two-point function is needed to 
compute ([O]), since higher-point correlators needed in principle have been expressed in terms 
of F{x,kj^), the cross section is still a non-linear function of that gluon distribution, invalidating 
-factorization. The rather simple form of the ^7 -factorization breaking term is due to the use 
of a Gaussian CGC color source distribution, and to the \wge.-Nc limit. However, the validity of 
this approximation has been critically examined |^], as it does not allow to correctly implement 
the non-linear QCD evolution, even in the large A'^,. limit. Finally, let us comment on the factor 



5(xP+— in formula (1.4). This delta function is a manifestation of the fact that in a high- 
energy hadronic collision, the momentum transfer is mainly transverse, and it appears because 
the eikonal approximation was used to compute the qA — qgX cross section. This is valid in the 
high-energy limit, as for instance the energy loss of the incoming quark is neglected. 

The CGC is endowed with a set of non-linear evolution equations which in the \wge.-Nc limit 
reduce to the Balitsky-Kovchegov (BK) equation [|l0|]. These equations can be interpreted as a 
renormalization group equation for the x evolution of the unintegrated gluon distribution, and more 
generally of ?i-point correlators, in which both linear radiative processes and non-linear recombi- 
nation effects are included. In this work, we compute the small-.x dynamics of the dipole correlator 
^ by solving the running-coupling (rc) BK equation. The evolution kernel is evaluated accord- 
ing to the prescription of Balitsky. Explicit expressions, together with a detailed discussion on 
the numerical method used to solve the rcBK equation can be found in JTl]], along with detailed 
discussions about other prescriptions proposed to define the running-coupling kernel. 



The only piece of information left to fully complete all the ingredients in ( |1.2| ) are the initial 
conditions for the rcBK evolution of .yV{x,r). This non-perturbative input has been constrained 
by single-inclusive forward hadron production data in [^]. The two parameters are xq = 0.02, the 
value of Xa below which one starts to trust, and therefore use, the CGC framework, and the value 
of the saturation scale at the starting point of the evolution 2^(xo) = Q^^ = 0.4 GeV^. Then, this 
information was simply taken over in [|5|]. In this respect, these forward di-hadron calculations aie 
predictions, there are no free parameter to play with. 
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Figure 1: The coincidence probability at a function of A^. Left: CGC calculations [||] for p+p and central 
d+Au collisions, the disappearance of the away-side peak is quantitatively consistent with the STAR data. 
Right: CGC predictions for different centralities of the d+Au collisions, the near-side peak is independent 
of the centrality, while the away-side peak reappears as collisions are more and more peripheral. 



2. Comparison with data 

We will now investigate the process dAu^h\h2X, with ^/s = 200 GeV. In particular we shall 
study the A0 dependence of the spectrum, where is the difference between the azimuthal angles 
of the measured forward particles hi and h2. To be more specific, we shall compute the coincidence 
probability to, given a trigger particle in a certain momentum range, produce an associated particle 
in another momentum range. It is given by 



CP{A(p) 



Npair{A(l>) 



Nt 



with NpairiA(t)) 



trig 



d^p\d^P2 



N, 



trig ■ 



yi^Pi±\ 



(2.1) 



In order to compare with the STAR measurement, the integration bounds for the rapidities are set 
to 2.4 < J < 4, which also ensures that only small-momentum partons are relevant in the nucleus 
wave function. In addition, for the trigger (leading) particle > 2 GeV and for the associated 
(sub-leading) hadron 1 GeV < \p2i_\ < \p\l\ - The single-inclusive hadron production spectrum, 
used to normalize the coincidence probability, is calculated as explained in [^. 

To deal with the centrality dependence, we identify the centrality averaged initial saturation 
scale 2^Q, extracted from minimum-bias single-inclusive hadron production data, with the value of 



i2 .^t u — 5 47 fjjj^ j-ijg Woods-Saxon distribution 7a(^) to calculate the saturation scale at 

other centralities: 

.2,.^_ QlMb) ^2 _o.4GeV^ 



(2.2) 



7^(5.47 fm) ' 

The result for central d+Au collisions is displayed in Fig.l, left plot, along with preliminary data 
from the STAR collaboration. As mentioned before, we do not calculated the complete near-side 
peak, as our formula does not apply around A0 = 0. We see that the disappearance of the away-side 
peak in central d+Au collisions, compared to p+p collisions, is quantitatively consistent with the 
CGC calculations. The latter are only robust for central d+Au collisions, but the extrapolation to 
p+p collisions is displayed in order to show that it is qualitatively consistent with the presence of 
the away-side peak, and also with the fact that the near-side peak is identical in the two cases and is 
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not sensitive to saturation physics. Since uncorrelated background has not been extracted from the 
data, the overall normalization of the data points has been adjusted by subtracting a constant shift. 

In Fig. 1, right plot, we show the centrality dependence of the coincidence probability. Al- 
though it is difficult to trust our formalism all the way to peripheral collisions, we predict that the 
near-side peak does not change with centrality, and that the away-side peak reappears for less cen- 
tral collisions. This is consistent with the fact that peripheral d-i-Au collisions are pH-p collisions. 
The fact that the away-side peak disappears from peripheral to central collisions shows that indeed 
this effect is correlated with the nuclear density. Moreover di-hadron correlations at mid-rapidity, 
which are sensitive to larger values of xa, feature an away-side peak whatever the centrality. The 
fact that for central collisions the away-side peak disappears from central to forward rapidities also 
shows that the effect is correlated with the nuclear gluon density. In a similar way, we predict 
that for higher transverse momenta, the away-side peak will reappear, as larger values of xa will 
be probed. These modifications of the away-side peak with transverse momenta and rapidity were 
already predicted in [|8|] as a time where there was no data. 

We are not aware of any descriptions of this phenomena that does not invoke saturation effects. 
We note that apart from our CGC calculation, a successful description based on the KLN saturation 
model was also recently proposed [|^. While more differential measurements of the coincidence 
probability, as a function of transverse momentum or rapidity, will provide further tests of our CGC 
predictions and help understand better this theory of saturation, the analysis of forward di-hadron 
correlations presented in this work adds further support to the idea that the saturation regime of 
QCD has been probed at RHIC. Future p-i-Pb collisions at the LHC will allow definitive tests. 
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